arXiv: 1503.04234v2 [nlin.PS] 25 Sep 2015 


Coexistence of stable branched patterns in anisotropic inhomogeneous systems 


B. Kaoui 1 , A. Guckenberger 1 , A. Krekhov 1,2 , F. Ziebert 1,3 and W. Zimmermann 11 

11 Theoretische Physik, Universitat Bayreuth, 95440 Bayreuth, Germany 
2 Max-Planck-Institute for Dynamics and Self-Organization, 37077 Gottingen, Germany 
3 Physikalisches Institut, Albert-Ludwigs-Universitat Freiburg, 79104 Freiburg, Germany [j 

(Dated: September 28, 2015) 

A new class of pattern forming systems is identified and investigated: anisotropic systems that 
are spatially inhomogeneous along the direction perpendicular to the preferred one. By studying the 
generic amplitude equation of this new class and a model equation, we show that branched stripe 
patterns emerge, which for a given parameter set are stable within a band of different wavenum¬ 
bers and different numbers of branching points (defects). Moreover, the branched patterns and 
unbranched ones (defect-free stripes) coexist over a finite parameter range. We propose two systems 
where this generic scenario can be found experimentally, surface wrinkling on elastic substrates and 
electroconvection in nematic liquid crystals, and relate them to the findings from the amplitude 
equation. 


I. INTRODUCTION 

Pattern formation is one of the most fascinating and intriguing phenomena in nature [HE]. It takes place in a 
wide variety of physical, chemical and biological systems and on disparate spatial and temporal scales, for example, 
convection phenomena in geoscience (HE] or in liquid crystals m, environmental patterns on, or patterns in 
chemical reactions uni uu and bacterial colonies [12 . In some circumstances pattern formation is undesired, for 
instance, the formation of spiral waves leading to cardiac arrhytmias in the heart muscle m- In other contexts 
pattern formation is even essential for the functioning of a system, e.g. in embryo development m or when designing 
surface wrinkling patterns to fabricate nanometer-scale structures pa ng. The mechanisms leading to the same 
type of pattern in different systems are obviously very diverse. Nevertheless, patterns occuring in systems of the 
same symmetry share common qualitative properties that can be described by universal amplitude equations for the 
envelope of periodic patterns EJISIISIEEZ]. 

Generating, modifying or eliminating patterns hence either requires a profound understanding of the pattern forma¬ 
tion mechanism in each specific system, or complementary, of the universal properties of a class of patterns and their 
response to symmetry breakings. Common pattern interventions are feed back control m and symmetry breaking via 
spatial, temporal or spatio-temporal modulations [T91I3Q] . To mention two interesting scenarios, spatial forcing near 
resonance (between the forcing and the natural wavelength) can lead to so-called incommensurate patterns [19j 20], 
while symmetry breaking via long-wave spatial modulations can render stationary patterns time dependent pTH34j . 
The response of patterns in quasi one-dimensional (ID) systems is now well established, however truely 2D scenarios, 
like the interplay of an anisotropy and a modulation in different directions, are yet fully unexplored. 

In this paper we identify and analyze a new class of quasi-2D anisotropic inhomogeneous pattern forming systems. 
The wave vector of the patterns lies close to q 0 = (go,0) along the preferred x-direction (anisotropy), as for the two 
experimental systems sketched in figure [l] (a) a spatially modulated version of the wrinkle forming system unnssi and 
(b) a modulated version of electroconvection (EHC) in nematic liquid crystals [6]. By a modulation we break in such 
systems the translational symmetry along the perpendicular ^/-direction, causing a variation of the pattern’s natural 
wavenumber go- This can be accomplished by varying the elasticity in the wrinkle forming system or the height of the 
electroconvection cell, respectively. In this class of anisotropic systems we find straight stripes, cf. figure [lj(c), that 
are stable for wavenumbers in a finite range around go, similar as in homogeneous anisotropic systems. In addition, 
however, we find a whole family of stable branched patterns as shown in figure [ljd) and (e). Surprisingly, they 
have - at identical parameters - different characteristic wave numbers and include different numbers of branching 
points. Moreover, the branched patterns coexist with the straight stripes in a wide parameter range. This behavior 
is a non-trivial generalization of the wave number bands (Eckhaus bands) for homogeneous systems [371443] - a well 
established concept and experimentally verified e.g. in EHC [35] and axisymmetric Taylor vortex flow 0H [42] - to 
inhomogeneous systems and multiple patterns. In the following we present the universal amplitude equation of this 
new symmetry class of patterns and analyze its solutions. 
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FIG. 1: In (a) a wrinkle forming elastic system is sketched with a spatially varying stiffness and in (b) electroconvection in 
nematic liquid crystals with a spatially varying height of the convection cell. Parts (c)-(e) show selected characteristic patterns 
for the proposed system class. They are stable at identical parameters and have been obtained numerically via equation § 
for the parameter set go = 1, k m = 0.05, M = 0.1, e — 0.08 and for different initial conditions by using the pseudo-spectral 
method [35]. Each picture is a cutout of a larger domain (L x — 80 x 27r/go and L y — 2n/km)- 


II. MODEL AND GENERIC AMPLITUDE EQUATION 


A generic model for the formation of stationary periodic patterns in anisotropic - but homogeneous - 2D systems, 
described by a field u(x, y , £), has been proposed in Ref. [Hj . One interpretation of the field u(x, y , t) is, that it describes 
the (small) lateral displacement of a thin elastic plate extended in the x-y plane, loaded along the x and y direction 
and supported by an elastic medium, similar as in wrinkling systems (whereby the in-plane elastic deformations are 
neglected) [43 . Here we generalize this model to an inhomogeneous situation by modulating the pattern’s preferred 
natural wavenumber qo along the direction perpendicular to the anisotropy (the x-direction), 

Qo(y) = % + Mcos{k m y), (1) 


with an amplitude M and a wavenumber k m considerably smaller than go- Then the dynamics of the patterns, 
described by the real field u(x,y,t) is goverend by 


d t u = 


ko(y) + A ) 2 


u — Wd^dyii — cdyU — u 3 


- 2{dyu)d v [<?o(2/)] -udy[qo(v)\ ■ 


( 2 ) 


The first line corresponds to the original model in Ref. (43] and the second line includes contributions due to the 
modulation qo{y). Equation is a representative of the here-identified symmetry class. It can be directly linked 
to the elastic wrinkle-forming system via an appropriate rescaling [43]: q$ = relates the critical wave number 
to the bending stiffness Ai of the hard layer and the elastic modulus k of the substrate. The control parameter 

£ - if) ”» is related *° * he critical MM = Abo note that, similar to the homogeneous 
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version [43] and the wrinkle system, equation (J2| can be derived from a functional as described in the appendix. 
For the following, we have chosen W = 1, c = 0.5 and go = 1 that favor straight stripes formation parallel to the 
^/-direction for the homogeneous case. 

Close to the onset of supercritically bifurcating patterns their generic (system-independent) properties may be 
described by a nonlinear dynamical equation for the complex envelope A(x,y,t), considering u(x,y,t) = Ae l ^° x + 
A* e ~iQo x w ith the scaling u oc e 1 ^ 2 (the star denotes the complex conjugate). This reduction method, the so-called 
multiple scale analysis, is well established for supercritical bifurcations in 2D homogeneous isotropic mm and 
anisotropic systems um or ID inhomogeneous ones [32]. Following very closely the appendix in Ref. [2] with the 
same intermediate scaling for time and space one obtains with M oc s 1 / 2 and k m oc s 1 / 2 via the multiple scale analysis 
from the system described by equation © the generic amplitude equation for the here-identified universality class of 
patterns in anisotropic systems with a spatially varying natural wave number qo(y)' 

d t A =[e + £dl + £ 2 d 2 ] A - g\A\ 2 A 

- *8<?qM cos(k m y)d x A - 4q^M 2 cos 2 (k m y)A. (3) 

Upon derivation from equation © one obtains the relations £ 2 = 4q^ , £ 2 = Wq^ and g = 3. The generic phenomenon 
of stripes’ branching, as detailed below, is induced in equation © by the term oc Md x A. The term oc M 2 A mainly 
leads to quantitative modifications. It should be noted that the first line of equation © is exactly the equation for the 
amplitude of stripe patterns in unmodulated EHC mm , the second anisotropic system suggested for experimental 
investigations of phenomena described here. For all specific pattern forming systems in the considered universality 
class, equation © can be derived via perturbation techniques from the basic equations (as for example from equation 
(J2|) and it has always the same form — only the coefficients will depend on the specific system. The presented results 
are obtained via the universal equation ©• 


III. RESULTS AND DISCUSSION 

The modulation of the natural wavenumber, see equation ([!]), alters the bifurcation from the basic state (u = A = 0) 
towards stationary periodic solutions of the form A(x,y) = e l ® x H(y). The stability of the latter and of A = 0 
with respect to small perturbations v(x,y,t) = vi(y)e at + lKx + (y)e a t ~ lKx can be determined via the ansatz 
A(x,y,t) = z lC * x [H(y) + v(x,y,t)\, followed by a linearization of equation ^ with respect to small v. An expansion 
of H(y) and (y) in Fourier modes results in an eigenvalue problem for cr, whereby solutions are stable when all 
eigenvalues have negative real parts. Note that in v(x,y,t) we consider only perturbations along the x-direction, as 
perturbations parallel to the stripes (along y) have been numerically found to have no effect. 




FIG. 2: (a) For small modulation amplitude, M = 0.03, the neutral curve (solid), above which stripes exist, and the Eckhaus 

stability region (in red), where stripes are stable, have the classical single-well shape with a minimum at Q = 0. For M — 0 
the neutral (dashed) and the Eckhaus (dotted) curves are given for comparison, (b) With M — 0.075 beyond a critical value 
(here M c ~ 0.0459), the neutral curve develops a double-well shape (solid line). The Eckhaus stable region for straight stripe 
patterns (in red) shrinks and detaches from the neutral curve: a gap develops, where only branched stripes are stable. Within 
the red region, stable unbranched and branched stripes coexist, (c) The minima of the neutral curve, Qmin , and (d) the gap 
8e between the neutral and the stable range as a function of M. Other parameters: go = 1 and km =0.1. 
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FIG. 3: (a) The amplitude H(y) of a modulated stripe pattern given in (b). (c) The energy density difference between the 

stripe pattern and one of the possible branched patterns, 5f = f s — fa, as a function of e for three different wavenumbers [Q = 0 
(solid), Q = 0.0125 (dashed), Q = 0.025 (dotted)] within the coexistence region. For small values of e, branched patterns are 
energetically preferred, while for large e unbranched ones. The crossover position strongly increases with Q. Parameters: 
qo = 1, km = 0.1, M — 0.075. 

For the basic state A = 0 the stability condition, 5ft [o- max (£ : Q)\ = 0, determines the neutral curve £n(Q), above 
which stationary periodic solutions, A = e l Q x H(y), of finite amplitude do exist. The neutral curves for the unmod¬ 
ulated case M m 0 and for a small modulation M = 0.03 are plotted in figure [2ja) with dashed and solid lines, 
respectively. Both curves have their minimum at Q = 0 and are symmetric. However, when increasing M, the value 
of the minimum £n(0) shifts towards positive values due to the term oc M 2 A in equation ^ that becomes larger 
(e —)> [s — Aq^M 2 cos 2 (k m y )]). By increasing M further beyond a critical modulation (here M c = 0.0459 for % = 1 
and k m = 0.1) a very important phenomenon takes place: The neutral curve develops two minima at finite wave 
numbers Q = ±Q m ^ n , as shown in figure [2|b) by the solid line. The appearence of these two minima is caused by 
high modulation amplitudes and mainly by the term oc Md x A in equation The dependence of Qmi n oc y/M — M c 
on the modulation in the vicinity of M c resembles a pitchfork bifurcation, cf. figure [2|c). The two new terms in 
the amplitude equation, caused by the modulation, hence widen the neutral curve at small modulations and change 
its shape from a single- to double-well shape via a pitchfork-like bifurcation at higher modulations. The fact that 
the modulation along the y-direction causes two minima of the neutral curve £n(Q)i with respect to Q along the x 
direction, is crucial for the emergence of branched patterns. 

Periodic solutions in the homogeneous case, M = 0, have a constant amplitude, H(y) = y/(e — i 2 Q 2 )/g , and they 
are linearly stable only above the dotted line in figure [2^a), the so-called Eckhaus-stability boundary £e(Q )• Note 
that for 2D anisotropic systems, the Q-range of stable solutions is symmetric with respect to Q = 0 [43] . in contrast 
to isotropic systems, in which the zig-zag instability occurs for Q < 0. For small modulations M < M c , H(y) is 
27r/£; m -periodic along the ^/-direction and corresponds to modulated stripes as shown in figures [lja) and[3ja ),(b). In 
general, the equation for H(y) may have further solutions as well. However, simulations did not show them, so they 
are either unstable or have higher energies. The Eckhaus stability range of the modulated stripes is qualitatively 
unchanged compared to the homogeneous limit, as shown by the red-colored region in figure |2ja): It still touches 
£n{Q ) at Q = 0. However, quantitatively, the modulation tends to narrow the width of the Eckhaus-stability range 
and thus reduces the range of stable stripes. 

Concomitant with increasing the modulation beyond the threshold M c a striking change in the stability scenario 
occurs: i) As shown in figure |2|b), a gap opens between the neutral curve (solid line) and the stability range of 
amplitude modulated straight stripes (red-colored area). Withhin the gap, the stripe solutions A(x,y) = e l ® x H(y) 
are unstable and thus cannot be observed anymore. The width of this gap increases nonlinearly with M — M c , as 
shown in figure [2jd) . ii) Stable branched stripes, like those shown in figures CTd),(e) and figure |4j emerge within a 
large area above the neutral curve, iii) Within the red-colored area in figure [2pc>), the unbranched modulated stripes 
coexist with branched stripes at identical parameters sets. 

To characterize the coexistence (the first central result) further, we use the fact that the dynamical equations © 
and © can be derived from functionals, which for stationary patterns take the simple form (appendix): 

T — — - J dxdy u 4 = — ^ J dxdy \A\ 4 . (4) 

Using the energy density / = T/L x L y (L x L y is the domain area) one can hence determine which of the stable 
patterns is energetically preferred at a given parameter set (see also p. 868 of Ref. El)- We have calculated the 
energy densities of unbranched patterns f s and of branched patterns fa within the red-colored region in figure [2jb) . 
Their difference Sf = f s — fa is plotted in figure [3^c) as a function of 6 and for three different wave numbers Q. 
For small 5, branched patterns have a lower energy and thus are energetically preferred, whereas for considerably 
larger 5, unbranched patterns are in turn preferred. The crossover can be understood as the modulation becomes less 
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FIG. 4: The branched patterns display different arrangements of defects upon increasing the amplitude M: (a) Aligned 
structure for M — 0.06, (b) Zig-zag order for M — 0.1, (c) Cascade-like order for M = 0.3. Parameters: qo = 1, km = 0.05, 
e — 0.08, L x — 80 x 27r/^o, L y — 2jr/km- 

important at large values of the control parameter 6. 

The second central result is that for each parameter set out of a large range beyond £n{Q ) we find always a whole 
family of stable branched patterns, characterized by different wavelengths and distinct numbers of branching points - 
and obtained by varying the initial conditions. For example, the two patterns in parts (d) and (e) in figure [l] contain 
five and seven pairs of branching points at identical parameters. They are both stable, however, the value of the 
functional of every branched pattern is, in general, different. This pattern coexistence in inhomogeneous systems is a 
surprising generalization of the well established Eckhaus wave number band of stable stripe patterns in homogeneous 
systems [37H42] . There, stable stripes of different wave numbers coexist at an identical parameter set, as verified in 
EHC [39], buckling [40] and axisymmetric Taylor vortex flow mm- Here, a band of branched patterns with different 
wavelengths and numbers of branching points coexists. In both cases, the origin of this coexistence lies in the fact 
that transitions between two stable periodic patterns require rather high excitations [38] worth to be investigated 
further. 

Another interesting observation is due to the branching points (defects) having opposite topological charges in each 
half period tt/ k m along the ^/-direction hence they repel each other. As shown in figure [IJ this repulsion can lead to 
different orderings of the defects with increasing defect density, for instance, the zig-zag ordering shown in figure |4][b). 
It can even lead to more complex cascade-like ordering, as shown in figure |4][c). However, how the detailed ordering 
of branching points can be controlled by the magnitude, the wavelength and the anharmonicity of the modulation is 
an interesting and fundamental question that needs to be adressed in the future. 

The phenomena identified here are robust and insensitive to the special shape of the imposed modulation of 
the natural pattern wavenumber. Studying the generic model, equation and imposing instead of a harmonic 
modulation either an anharmonic wave-number modulation or a step-wise one, we obtain qualitatively the same 
results [35]. Consequently, instead of harmonic variations, which may be challenging to implement in experimental 
systems, a step-like modulations can be used, which is easy to achieve in the wrinkles system by glueing two elastomer 
substrates with different elastic properties together m m- 


IV. CONCLUSIONS 

In conclusion, we identified a new symmetry class of pattern forming systems that are anisotropic with a perpen¬ 
dicular modulation. By deriving and analyzing the universal amplitude equation, we showed that this class displays 
interesting scenarios - the emergence of families of stable branched patterns and their coexistence with unbranched 
patterns - that we suggest to verify experimentally in two complementary anisotropic systems: First, in wrinkle 
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forming systems m, where the elasticity of the substrate [as sketched in figure [lja)] and/or the thickness of the hard 
layer on top can be varied perpendicularly to the pre-stretch (anisotropy) direction, as recently realized experimen¬ 
tally [35, 44. And second, in dissipative electroconvection in nematic liquid crystals mg, where the layer height 
[as sketched in figure [ljb)] or the driving frequency are modulated perpendicularly to the nematic alignment. For 
both systems, a thorough theoretical analysis of the modulated basic equations is feasible, including the derivation 
of the universal equation of the patterns envelope, i.e. equation ©>• There is a vast literature on wrinkles recently 
(see e. g. Refs. pH M EH HIM]), including control strategies for wrinkle formation, but the question of pattern 
coexistence seems to have not even been raised yet. Hence studies along the lines proposed here, including quantitative 
comparisons between theory and experiments, will prove powerful for future control and design strategies, both of 
unbranched and branched (wrinkle) patterns. Finally, note that the branching of stripe patterns on the skin of fishes 
[48] . an important example for patterns in living systems, is probably also induced by inhomogeneities, similar as 
discussed here, and which are in addition slowly changing with time. 
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Appendix A: Functionals 


The dynamical equation Eq. © for the field u(x,y,t) can be derived from the functional 


F[u{x,y,t)\ =J f{u)dxdy= ^ J 


-eu 2 + ^ u 4 + W ( d x d y u ) 2 


via the functional derivative 


, , >F 8/ , 3/ 

d t u(x,y,t) = ~ r -\ Yu +d.. 


+c (d 2 u ) 2 + (qo(y)u 2 4- 2q%(y)uAu + (A u) 2 ) 

df 


dx dy 


: d(d%u) 


d x dy —t + d 2 ^ 


d(d x d y u ) y d(d?u) / ' 


(Al) 


(A2) 


The functional in Eq. © can be simplified by assuming periodic boundary conditions and by using the expression 
d t u from the dynamical equation. One gets 


F[u(x,y,t)\ =4/ 


2 u4 ~ u ( d tu) -u 4 - 2 u(d y u)dy{ql(y)) - u 2 d 2 [q 2 (v)] 


dx dy , 


(A3) 


where the — rU-term compensates for the cubic term occurring in d t u and similarly the last two terms. These last two 
terms again vanish due to the periodic boundary condition (here in ^/-direction): 


J [~ 2u (9yu)d y (qo(y)) - u 2 &l [<7o(y)]] dxdy = -Jd y [u 2 d y (< ql{y ))] dxdy = 0. 
We arrive at the simple expression 


(A4) 


^u 4 + u(d t u) 


dx dy 


T[u(x,y,t)\ = ~\ J 

and if we are - as in our paper - interested only in stationary patterns, we simply have 

^[^(x, 7/, t)] = — - J u A dxdy . 

Also the amplitude equation |3| can be derived from a functional, which can be simplified for stationary solutions ir 
the same manner, only the prefactor changes: 


(A5) 


(A6) 


T[u{x,y,t)} = - ^ J \ A\ 4 dxdy. 


(A7) 
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